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ABSTRACT 

We report on the results from observations of the most recent outbursts of XTE J1739— 302 
and IGR J17544— 2619, which are considered to be the prototypes of the supergiant fast X-ray 
transient (SFXT) class. They triggered the Swift/BAT on 2011 February 22 and March 24, 
respectively, and each time a prompt Swift slew allowed us to obtain the rich broad-band 
data we present. The XRT light curves show the descending portion of very bright flares that 
reached luminosities of ~ 2 x 10 36 and ~ 5 x 10 36 erg s _1 , respectively. The broad-band 
spectra, when fit with the usual phenomenological models adopted for accreting neutron stars, 
yield values of both high energy cut-off and e-folding energy consistent with those obtained 
from previously reported outbursts from these sources. In the context of more physical models, 
the spectra of both sources can be well fitted either with a two-blackbody model, or with a 
single unsaturated Comptonization model. In the latter case, the model can be either a classical 
static Comptonization model, such as COMPTT, or the recently developed COMPMAG model, 
which includes thermal and bulk Comptonization for cylindrical accretion onto a magnetized 
neutron star. We discuss the possible accretion scenarios derived by the different models, and 
we also emphasize the fact that the electron density derived from the Comptonization models, 
in the regions where the X-ray spectrum presumably forms, is lower than that estimated using 
the continuity equation at the magnetospheric radius and the source X-ray luminosity, and we 
give some possible explanations. 
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X-rays: individual: XTE J1739-302, IGR J17544-2619. 



1 INTRODUCTION 

Supergiant fast X-ray transients (SFXTs) are a class of High Mass 
X-ray Binaries (HMXBs) associated with OB supergiant stars. In 
the X-rays they display outbursts significantly shorter than those 
of typical Be/X-ray binaries characterized by bright flares with 
peak luminosities of 10 36 -1 37 erg s -1 which last a few hours 
(as obs erved by INTEGRAL; ISguera et al]|2005l : iNegueruela et al.l 
2006b). As their qui escence is characterized by a lum inosity of 
- 10 32 erg s" 1 (e.g. lin't Zanj|20ol : lBozzo etaifeoioh . their dy- 
namic range is of 3-5 orders of magnitude. While in outburst, their 
hard X-ray spectra resemble those of HMXBs hosting an accreting 
neutron star (NS), with hard power laws below lOkeV with high 
energy cut-offs at ~ 15-30 keV. So, even if pulse periods have 
only been measured for a few SFXTs, it is tempting to assume that 
all SFXTs might host a NS. 



The physical context originating the outbursts has been 
claimed to be related eith er to the properties of the wind from 
the supergiant companion (j in't Zand 2005; Walt er & Zurita Herasl 
Negueruela et a l. 2008; Si doli et alj|2007l) or to the presence 



the sup* 
l2007l; fr 



of a centrifugal o r magnetic barrier (IGrebenev & Sunvaevll2007l: 
iBozzo et al]|2008h . 



XTE J1739- 302 was discovered by RXTE in August 1997 
dSmithetal.lll998h . when it reached a peak flux of 3.6 x 10 erg 
cm -2 s _1 (2-25 keV). It has a lo n g history of flaring recorded 
by INTEGRAL dS guera et al.l l2006t Iwalter & Zurita Herasl 120071; 
i Blav et al.ll2008h an d by Swift dSidoli et al.ll2009bl lal; lRomano et all 
uOllcl lbr). Recently. iDrave et al.l (2010) reported the discovery of 
a 51.47 ± 0.02 d orbital period based on ~ 12.4Ms of INTE- 
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Table 1. Log of the observations of XTE J1739-302 and IGR J17544-2619 with Swift. 



Source Sequence Instrument Start time (UT) End time (UT) Exposure Time Since Trigger 



/Mode (yyyy-mm-dd hh:mm:ss) (yyyy-mm-dd hh:mm:ss) (s) (s) 



XTEJ1739-302 00446475000 BAT/evt 2011-02-22 07:19:44 2011-02-22 07:36:46 1022 -119 

00446475000 XRTAVT 2011-02-22 07:24:10 2011-02-22 09:04:37 1271 147 

00446475000 XRT/PC 2011-02-22 09:04:39 2011-02-22 09:13:55 554 6176 

IGR J17544— 2619 00449907000 BAT/evt 2011-03-24 01:53:04 2011-03-24 02:13:06 1202 -239 

00449907000 XRTAVT 2011-03-24 01:59:15 2011-03-24 03:18:22 657 133 

00449907000 XRT/PC 2011-03-24 02:10:07 2011-03-24 04:01:03 3360 784 

00035056150 XRT/PC 2011-03-27 19:54:23 2011-03-27 20:11:56 1046 323840 

00035056151 XRT/PC 2011-03-28 20:10:36 2011-03-28 20:21:58 662 411213 

00035056152 XRT/PC 2011-03-29 20:15:13 2011-03-29 20:23:57 504 497890 



GRAL datfl The optica l counterpart is an 081 ab star at 2.7 kpc 
dNegueruela et al.ll2006al;lRahoui et alj|2008l) . 

IGR JT7544-2 619 was first detected by INTEGRAL in 2003 
dSunvaev et al.l2003h . when the source reached a flux of 160 mCrab 
(18-25 keV). Several more flares, lasting up to 10 hours, were 
detec te d by I NTEGRAL in the fo l lowing years dGrebenev et all 



20031 12004 



Kuulkers et al 



Sguera et alj l2006t IWalter & Zurita Herasl 120071: 
20071) with fluxes up to 400 mCrab (20-40 
keV); some were al s o fou nd in archival BeppoSAX observa- 
tions din't Zand et al, I 120041) . Subsequent flares were observed 
by Swift |Krimm et alj 12001 ; ISidoli et alJl2009bl P: IRomano et all 

"iJr 



I201ldl3lal) . and Suzaku dRampv et alj2009l) . IClark et alj d2009h re- 
ported the discovery of a 4.926 ± 0.001 d orbital period base d on 
~ 4.5 years of INTEGRAL data. Recently. iDrave et alj $20V%) de- 
tected a transient 71.49 ± 0.02 s signal in RXTE observations from 
the region of IGR J 17544— 2619, which, if interpreted as the spin 
period of the NS in the system , places the source in the locus of the 
Corbet diagram dCorbej(l986h where classical wind-fed supergiant 
X-ray binaries can be found. The optical counterpar t is an 091b 
star at 3.6 kpc dPellizza et al.ll2006tlRahoui et alj|2008l) . 

XTE J1739-302 and IGR J17544— 2619 are considered the 
prototypes of the SFXT class, and were extensively studied with 
Swift. In particular, in addition to the Burst Alert Telescope (BAT, 
iBarthelmv etai] l2005h outburst dete ctions and intensive X-ray 
Telescope (XRT, Burrows et alj|2005 ) follow-up, Swift has been 
studying their long-term properties dRomano et alj|2011cl and ref- 
erences therein). In this paper we examine the most recent outbursts 
of these two sources, which triggered the BAT in 2011. In partic- 
ular, we apply, for the first tim e to SFXTs, the new COMPMAG 
model bv lFarinelli et alj d2012l hereafter F12). 



2 OBSERVATIONS AND DATA REDUCTION 

XTE J 1739 -302 triggered the BAT on 2011 February 22 at 
07:21:37 UT (image trigger=446475, IRomano et al1l2011bl) . Swift 
immediately slewed to the target, so that the narrow-field instru- 
ments (NFI) started observing about 141 s after the trigger. After 
the initial automated target (AT) observation (two orbits for a to- 
tal of ~ 1.8 ks net exposure spanning ~ 6.7 ks) no further NFI 
observations were performed. 



1 We note, however, that IRomano et al ] 120091) derived a marginal evi- 
dence, based on XRT data of signal a t P olb = 12.8658 ± 0.0073 d, or 
1/4 the period reported bv lDrave et al.l J2010t) . 
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Figure 1. Swift/XKT (0.2-10 keV) long-term light curves of 
XTEJ1739-302 and IGR J17544— 2619. The downward-pointing ar- 
rows are 3cr upper limits. The upward-pointing arrows mark flares that 
triggered the BAT Transient Monitor on MJD 54414 and 54906. Data u p to 
MJD ~ 55261 (grey) were published in lRomano et alj d2009ll201 lddh . 



IGR J17544-2619 dRomano et al]|2011ah triggered the BAT 
on 2011 March 24 at 01:56:57 UT (image trigger=449907, 
IRomano et al.ll201 lal) . Swift immediately slewed to the target, so 
that the NFIs started observing about 126 s after the trigger. The 
AT (sequence 00449907000) ran for two orbits, until ~ 7.4 ks after 
the trigger). Follow-up target of opportunity (ToO) observations for 
a total of 2.2 ks were obtained (sequences 00035056150-152) after 
the source emerged from Moon constraint. The data cover the first 
6 d after the beginning of the outburst. Table[T]reports the log of the 
Swift observations used in this paper. 

The BAT data were analysed using the standard BAT anal- 
ysis software within FTOOLS in the HEASOFT package (v.6.11). 
Mask-tagged BAT light curves were created in the standard energy 
bands, and rebinned to fulfil at least one of the following condi- 
tions, achieving a signal-to-noise (S/N) of 5 or bin length of 10 s. 
Spectra were extracted from the whole event lists and within time 
intervals strictly simultaneous with the XRT. Response matrices 
were generated with BATDRMGEN using the latest spectral redis- 
tribution matrices. Survey data products, in the form of Detector 



© 2012 RAS, MNRAS 000,[TJ|JJJ 



SFXTs with Swift: outbursts of prototypes 3 



100 



1000 10 4 10 5 

Time (s since 01:56:57 UT) 



o 






(a) \ 


Oj _ 








rCoun 
0.1 


* 






x § 


r XTE J1739-302 




i 


100 1000 10 4 10 s 
Time (s since 07:21:37 UT) 


10" 


XRT Count Rate 
0.01 0.1 1 10 


r ICR J 17544-26 19 




(b)j 

• | 

t i 

1 ^ 1 1 1 i 



Figure 2. Swift/XRT (0.2-10 keV) light curves of the 2011 outbursts of 
XTEJ1739-302 and IGR J17544— 2619, as followed by Swift/XRT, re- 
ferred to their respective BAT triggers (2011 February 22, 2011 March 
24). Points denote detections, downward-pointing arrows 3cr upper limits. 
Where no data are plotted, no Swift data were collected. Vertical dashed 
lines mark time intervals equal to 1 day, up to a week. 
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Figure 3. XRT and BAT light curves of the initial orbit of data of the 201 1 
February 22 outburst of XTE J 1739— 302 in units of count s _1 and count 
s — 1 detector - 1 , respectively. The empty circles correspond to BAT in event 
mode, filled circles to BAT survey mode data. 



Table 2. Absorbed power-law spectral fits of XRT data of XTE J 1739 -302. 



Plane Histograms (DPH), are available, and were analysed with the 
standard BATSURVEY software. 

The XRT data were processed with standard procedures 
(XRTPIPELINE vO. 12.6), filtering and screening criteria by using 
FTOOLS. We considered both windowed-timing (WT) and photon- 
counting (PC) mode data, and sele cted event grades 0- 2 and 0-12, 
for WT and PC data, respectively feurrows et al.l 2005). When ap- 
propriate, we corrected for pile-up by determining the size of the af- 
fected core of the point s pread function (PSF) by comparing the ob- 
served and nominal PSF ( Vaugh an et al.l2006t) , and excluding from 
the analysis all the events that fell within that region. Background 
events were accumulated from source-free regions. For our timing 
analysis light curves were created for several values of signal-to- 
noise ratio (SNR) and number of counts per bin; all were corrected 
for PSF losses, vignetting and background. For our spectral anal- 
ysis, we extracted events in the same regions as those adopted for 
the light curve creation; the data were rebinned with a minimum 
of 20-50 counts per energy bin, as appropriate, to allow x 2 fit- 
ting, with the exc eption of low-count-statistics spectra, in which 
Cash JCashll 19791) statistics and spectra binned to 1 count per bin 
were used, instead. Ancillary response files were generated with 
XRTMKARF, to account for different extraction regions, vignetting, 
and PSF corrections. We used the latest spectral redistribution ma- 
trices in CALDB (20110915). For the luminosity calculation we 
adopted a distance of 2.7 kpc and 3.6 kpc for XTE J1 739— 302 and 
IGR J17544-2619, respectively l iRahoui et al1l2008l) . 

All quoted uncertainties are given at 90 % confidence level 
for one interesting parameter unless otherwise stated. The spectral 
indexes are parametrised as F v oc v~ a , where F v (erg cm -2 s _1 
Hz - ) is the flux density as a function of frequency u; we adopt 
r = a + 1 as the photon index, N(E) oc E~ T (ph cm -2 s~ x 
keV -1 ). 



Spectrum" AT* V Flux c L d x 2 /d°f 

C-stat (doff 



WT (147-1343) 2.35Io;i4 0.79+o og 1.8 1.8 681/563 
WT (4625-4697) 3.18+°;e5 L13 ^a27 13 lA $2/40 
PC (6176-6730) 6.44±*;g 1.42±{[H 0.02 0.03 73.4(83) 



a Spectrum (seconds since trigger). 
b In units of 10 22 cm -2 . 

c Observed 2—10 keV fluxes (10~ 9 erg cm -2 s — 1 ). 

d 2-10 keV luminosities in units of 10 36 erg s — 1 , at 2.7 kpc. 

e Cash statistics (C-stat). 

3 ANALYSIS AND RESULTS 

Fig. Q] shows the Swift/XRT 0.2-10 keV light curve of 
XTE 11739-302 and IGR 117544-2619 throug hout our 2007- 
2009 monitoring program (Ro mano et al.l |201 let and references 
therein) background-subtracted and corrected for pile-up, PSF 
losses, and vignetting. All data in one observation (l-2ks) were 
generally grouped as one point, except for outbursts, which show 
up as a vertical lines on the adopted scale. The data already 
presented elsewhere are in grey. The details of the 2011 outbursts 
of XTE 11739-302 and IGR 117544-2619, referred to their 
respective BAT triggers, are shown in Fig. [2] In both cases the XRT 
caught the descending portion of a very bright flare in great detail. 
There is a similarity between the two outbursts, and also with 
previou s outbursts of these sources [see figure 6 of iRomano et al.l 
fl201 la) , for a compendium of the best light cu rves collected before 
the currently presented ones and figure 9b-e of lSidoli et al1d2009al) 
for a comparison of different outbursts of the same source]. 

The details of the data selection and spectroscopic analysis on 
both sources are reported in the remainder of this Section. Here we 
summarise the models we adopted for the broad-band spectra. We 
first considered the phenomenological models typically used to de- 
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Table 3. Best-fit parameters of the phenomenological models for XTE J1739— 302 and IGR J17544— 2619 broad-band spectra, respectively. In the latter case 
an absorption edge was also included. CPL=cut-off powerlaw; HCP=powerlaw multiplied by exponential cut-off with e-folding factor. 



XTEJ1739-302 IGR J17544— 2619 



Parameter 


CPL 


HCP 


CPL 


HCP 


jV H (10 22 cm" 2 ) 


, Q7+0.17 


1 76+ ' 17 
1-' D -0.15 


n 7 ,+0.06 


1 m+ a06 

l- ui -0.05 


£ cdgc (keV) 






22.98+^2 




^cdgc 






,,J - 0.43 




r 


u - z;, -0.10 


42+ ' 09 
u -^-0.09 


04+ - 07 


n cq+0.02 
u - oy -0.04 


E c (keV) 




9 '56+ ' 99 




r fifi +0.54 
°- DD -0.51 


Ef (keV) 


Q in+0.95 
y - ly -0.80 


4 60+ ' 37 


7 50+ 54 

'■ JU -0.49 


14.32+5^ 2 


^2-10 keV 


2.4 


2.3 


3.1 


3.2 


^l-lOOkcV 


7.0 


7.1 


7.6 


9.0 


T h 

-^2-10 koV 


2.0 


2.0 


4.8 


5.0 


T b 
1-100 keV 


6.1 


6.1 


11.9 


14.1 


X 2 /dof 


364/243 


336/242 


344/296 


324/297 



a Average observed flux in units of 10 9 erg cm 2 s 
b In units of 10 36 erg s" 1 . 



scribe the X-ray emission from accreting pulsars in HMXBs, i.e., i) 
simple absorbed (WABS in XSPEC) power laws; ii) absorbed power 
laws with exponential cut-offs (CUTOFFPL, hereon CPL); Hi) ab- 
sorbed power laws with high-energy cut-offs (HIGHECUT, hereon 
HCP). The simple advantage of this approach is that these fits yield 
easy-to-compare estimates of fluxes and luminosities. The disad- 
vantage, however, is that little physical insight can be obtained from 
such fits, therefore we also considered the following physical mod- 
els: iv) absorbed generic Comptonization models (COMPTT) in dif- 
fusion approximation for a disc geometry without a dynamical bulk 
component, v) a combination of two blackbodies (BB) with differ- 
ent temperatures and radii (BBODY RAD+BBODYRAD ), sometimes 
used to fit the spectra of magnetars dlsrael et al.ll2008h ; vi) the new 
COMPMAG model, recently developed by F12. We refer the reader 
to that paper for a detailed description of the algorithm. Here we 
briefly remind the reader that COMPMAG is based on the solution 
of the radiative transfer equation for the case of cylindrical accre- 
tion onto the polar cap of a magnetized NS. The velocity field of 
the accreting matter can be increasing towards the NS surface, or 
it may be described by an approximate decelerating profile. In the 
former case, the free parameters are the terminal velocity at the NS 
surface, 0, and the index of the law /3(z) oc z~ v , while in the sec- 
ond case the law is given by /3(r) oc —t. The other free parameters 
of the model are the temperature of the BB seed photons, fcTbt,, the 
electron temperature and vertical optical depth of the Comptoniza- 
tion plasma kT c and r, respectively, and the radius of the accretion 
column ro, in units of the NS Schwarzschild radius. The different 
combinations of these parameters determine the steepness of the 
spectrum at high energies and the rollover energy position. 



However, given that the Comptonization spectra can be deter- 
mined by only three quantities (slope, cut-off position and normal- 
isation), for practical purposes it is generally not possible to keep 
free all the above mentioned parameters, a procedure which other- 
wise lead them to be completely unconstrained during the fitting 
procedure. In the following Sections we will report in more detail 
the procedure that we adopted in the spectral fitting. 



3.1 XTE J1739 302 

The XRT light curve (Fig.[2ji) shows the descending part of a flare, 
that started off at a peak exceeding 30 counts s~ then decreases 
to about 2.5 counts s~ . A second flare is observed in the second 
orbit, with a count rate in the range 10-17 counts s _1 . In the third 
orbit the count rate was ~ 0.2 counts s _1 . In Fig. [3] are reported 
the first orbit data of XTE J1739— 302 in several energy bands. The 
event-by-event mask-weighted light curves only show the initial 
flare, which started at < T — 70 s with a slow rise, well defined 
in both soft bands. The source was still bright when the BAT event 
data end, at T + 900 s. The BAT survey data cover the same time- 
span. 

We extracted an XRT spectrum for each orbit of data as de- 
tailed in Table [2] and we fit them with an absorbed power-law 
model. The results are summarized in Table [2] By linking the Nh 
across the spectra, we can detect Nh variations, while the photon 
index is consistent with being constant, as previously reported in 
ISidoli et aLl f2009bl lah. 

The BAT average spectrum (T — 119 to T+903 s) can be fit 
by a simple power-law model with a photon index F — 2.67~tg']^ 
(.Fis-sokcV = 2.8 x 10 -9 erg cm~ 2 s _1 ); the fit, however, is 
not good (x 2 /dof=51/37) and the residuals indicate that a curvature 
is present in the spectrum, which is best assessed in a broad-band 
analysis. 

In order to perform broad-band spectroscopy of the outburst, 
we extracted strictly simultaneous spectra from the XRT and BAT 
event lists, i.e., in the time interval 147-903 s since the BAT trig- 
ger (see Fig. [3}. We fit them in the 0.5-10 keV and 15-60 keV en- 
ergy bands for XRT and BAT, respectively. Factors were included 
in the fitting to allow for normalisation uncertainties between the 
two instruments, constrained within their usual ranges (0.9-1.1). A 
simple absorbed power-law model clearly yields an inadequate fit 
of the broad-band spectrum with a x 2 /dof=926/244. A significant 
improvement is obtained when considering the other phenomeno- 
logical curved models described above (see Table|3}. 

A more satisfactory fit is instead obtained with the model 
WABS x (BBODYRAD+BBODYRAD), whose results are reported in 
Table [4] while in Fig. [4] we show the deconvolved best-fit model 
and data. 

We also tested the model WABS x COMPTT, which yields sta- 
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Figure 4. Absorption-corrected unfolded EF(E) models and spectra, and 
residuals between the data and the model in units of cr for XTE J1739— 302. 
Upper panel: BBODYRAD+BBODYRAD. Lower panel: COMPMAG. 



Table 4. Best-fit parameters of the 
WABSX (BBODYRAD+BBODYRAD) for XTE J1739— 302 
IGR J17544— 2619 broad-band spectra, respectively. 



model 
and 



Parameter 


XTEJ1739-302 


IGR J17544-2619 


Nu (10 22 cm" 2 ) 


1 07+°- 09 


U - 5S -0.06 


kT£ h (keV) 


, oc+0.07 
1-»S_0.06 


9 3+o.07 

u -"-0.07 


Rbb.l (km) 


i ic+0.04 
1 ' 1J -0.05 


3 4 3+0.41 


fcT 2 b (keV) 


6 62+ ' 49 

aDZ -0.15 


3 4^+0.08 
-'■ HJ -0.07 


Rbb,2 (km) 


u+ l) - 01 


n 77 +0.04 
u -"-0.04 


pa 
2-10 keV 


2.2 


3.1 


^l-lOOkcV 


5.9 


7.7 


T h 

-^2-10 keV 


1.9 


4.8 


T h 

-^1-100 keV 


5.2 


12.1 


X 2 /dof 


285/241 


335/297 



a In units of 10 9 erg cm 2 s 1 . 
b In units of 10 36 erg s _1 . 



tistically similar results, and we report the main values of the model 
in Tabled 

Finally, we tested the COMPMAG model assuming a free-fall 
like velocity profile (77 = 0.5), two terminal velocities at the NS 
surface, = 0.05 and = 0.2, respectively, and NS albedo 
A = 1. The free parameters of the model during the fitting pro- 
cedure were the seed BB temperature fcTbt>, the electron temper- 
ature fcT c , and the vertical optical depth of the accretion column 
r. We note that the latter quantity in this case is a factor about 



Table 5. Best-fit parameters of the model WABSX COMPTT for 
XTE J1739-302 and IGR J17544-2619 broad-band spectra, respectively. 
As for the case of Table [3] an absorption edge was included in the model 
for IGRJ17544-2619. 



Parameter 


XTE J1739-302 


IGR J17544— 2619 


Nu (10 22 cm" 2 ) 


US1 -0.10 


n SR +0.07 


£cd gc (keV) 

''"edge 

fcT w (keV) 
kT c (keV) 
r 


1 ^4+007 
1 - JH -0.06 

8.83+iff 
361 +o:5i 

-0.79 


20.97+^ 2 
62+ 3 ' 5 
63+ - 07 

-0.16 

10 67 +o.7i 


pa 

^2-10 keV 


2.2 


3.0 


pa 

^l-lOOkcV 


6.8 


8.0 


u 2-\0 keV 


1.9 


4.7 


fb 

^1-100 keV 


5.9 


12.5 


X 2 /dof 


289/242 


323/295 



a In units of 10 9 erg cm 2 s 1 . 
b In units of 10 36 erg s _1 . 



Table 6. Best-fit parameters of the model WABSX COMPMAG for 
XTE J1739-302 and IGR J17544-2619 broad-band spectra, respectively. 
In both cases the fixed parameters are r) = 0.5, = 0.05, rn = 0.25, and 
A = 1. Also in this case, as in TablefS] an absorption edge was included in 
the model for IGR J17544-2619. 



Parameter 


XTEJ1739-302 


IGR J17544— 2619 


Nh (10 22 cm" 2 ) 


1 1X+ 012 
11S -0.09 


nfi8 +o.io 

U - 68 -0.O8 


£od gc (keV) 




21 71+ 1 - 84 

Z1 -'-0.95 


''"edge 




0.81_ 23 


kT bh (keV) 


1 63+ ' 10 
1 -°-'_0.17 


73+0.11 


kT c (keV) 


9-« as? 


3 47 +0.34 
•'■^'-0.02 


T 


n to+0.30 
u ""-0.03 


1 iq+0.09 
^^-O.IO 


pa 

^2-10 keV 


2.3 


3.0 


-Fai_l00 keV 


7.4 


8.1 


T h 

-^2-10 keV 


1.9 


4.7 


T h 

-^1-100 keV 


6.4 


12.7 


X 2 /dof 


292/241 


324/295 



a In units of 10 9 erg cm 2 s . 
b In units of 10 36 erg s -1 . 



1/1000 lower than classical optical depth because of the inclusion 
of an energy-independent correction term which reduces the Thom- 
son cross-secti on or for photons prop agating in the direction of the 
magnetic field lBecker & Wolfj d20oi hereafter BW07). 

We initially set the radius of the accretion column ro = 0.25, 
which corresponds to ~ 1 km for a NS star with mass M = 
1.4Mq; the results for the case = 0.05 are reported in Table 
[6] while in Fig.|7]we show the deconvolved best-fit model and data. 

No significant variation is observed in the \ 2 -value when in- 
creasing ro from 0.25 to 1. A similar statistical result is also ob- 
tained assuming = 0.2 (x 2 /dof=294/241). In this case, a slight 
decrease is observed both in the electron temperature and optical 
depth, although their values can be considered unchanged within 
errors (kT c = 6.5±% s 8 keV and r = 0.291°,;,", respectively). The 
electron temperature drop is not surprising, as in a thermal plus 
bulk Comptonization scenario the spectral slope and cut-off energy 
are dictated not only by the plasma temperature and optical depth, 
but also by the shape of the velocity field and maximum velocity of 
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Table 7. Absorbed 
IGR J17544-2619. 



power-law spectral fits of XRT data of 



Spectrum" 



Flux c 



x 2 /dof 
C-stat (dof) e 



WT (133-783) 
PC (784-1659) 
PC (4881-7418) 



1.00 
1.83 
2.33 



+0.06 
0.06 
+0.55 
0.49 
+0.47 
0.41 



0.69 
1.05 
1.39 



+0.04 
0.04 
+0.28 
0.27 
+0.25 
0.23 



3.0 

0.27 

0.06 



4.89 545/570 
0.46 30/33 
0.11 46/42 



a Spectrum (seconds since trigger). 
b In units of 10 22 cm -2 . 

c Observed 2-10keV fluxes (10 -9 erg cm -2 s" 1 ). 

d 2-10keV luminosities in units of 10 36 erg s — x , at 3.6kpc. 

e Cash statistics (C-stat). 



E- 

K 

X 



10 
1 

10 

1 

0.015 
0.01 
0.005 


0.01 
0.005 




J 1 1 1 1 1 1 1 1 1 1 1 1 1 

1 IGR J17544-2619 *. „ 
r 0.3-4 keV 


I I I I | I I I 5 


r 4-10 keV V - 


H — i — i — i — i — — i — i — i — = 




15-25 k'eV ' ' j 




T T f T ] "T T 1 T — 

25-50 keV j 



-500 







500 



1000 



1500 



Time (s since 01:56:57 UT) 



the accreting matter (see fig. 4 in F12). Thus, for given (observed) 
spectral slope and roll-over energy, the higher the values of 0, the 
lower the values of kT c , as the bulk Comptonization progressively 
increases its importance at expenses of thermal Comptonization. 
Even in this case, with = 0.2, the x 2 remains unchanged for 
varying ro from 0.25 up to 1. 

An important issue we want to emphasize is that for 
XTE J1739— 302 the value of r obtained with COMPMAG is less 
than 1 (see Tabled, and actually the higher the assumed value of 
the lower the value of r. 

As the COMPMAG model is based on the solution of the 
Fokker-Planck approximation of the radiative transfer equation, 
it generally holds for an optical depth r > 1. Our fits on 
XTE J 1739— 302 data yield values of r which are formally smaller 
than unity. In general terms, this would indicate that the diffusion 
along the direction of the column axis is not efficient, hence ei- 
ther the spatial diffusion of the photons is only efficient for photons 
propagating perpendicularly with respect to the column axis, or the 
system is significantly anisotropic due to the presence of the mag- 
netic field. In either case, the diffusion approximation would not 
be recommended, as inadequately describing the system (further 
details can be found in the Appendix). 

We note however that, in our case, an acceptable fit can be 
achieved with a value of the optical depth within a factor of 3 
of unity, so that the obtained physical parameters can be dis- 
cussed with reasonable accuracy. In particular, the electron temper- 
ature kT c results consistent within errors with that obtained with 
COMPTT and with the hotter BB temperature in the two-BB model 
(see Tables[4]and[5]l. 



3.2 IGR J17544 2619 

The XRT light curve of IGR J17544-2619 is reported in Fig. f2}), 
and shows the two bright flares observed in the first orbit; the first 
was caught (WT mode) in the descending part and reached a peak 
exceeding 50 counts s _1 ; the second (PC mode) started off at about 
T+750 s and reached ~ 5 counts s -1 . At about T+5000 s (second 
orbit) the source is still at ~ 1 counts s - . Fig.|5]shows the initial 
data of IGR J17544— 2619 in several energy bands. The BAT event- 
by-event mask-weighted light curves show a slow rise from T — 
200 s up to the first peak at T — s, followed by a few more flares, 
then follows the decaying shape of the XRT light curve. The BAT 
survey data cover a longer time-span, and show that the source was 
still detected in the soft band out to T + 1500 s. 



Figure 5. Same as Fig. fj] for the 2011 March 24 outburst of 
IGR J17544— 2619. The empty circles correspond to BAT in event mode, 
filled circles to BAT survey mode data. 



We extracted three XRT spectra at different times since the 
trigger and we fit them with an absorbed power-law model, as de- 
tailed in Table|7] The results are summarized in Table|7] By linking 
the Nh across the spectra, we can detect a harder- when-brighter 
trend and Nn variations, as previously observed by iRampy et al.l 
(2009). We note that this is the brightest outburst of this source 
recorded by Swift. 

The BAT average spectrum (T — 239 to T+963 s) was fit by 
a simple power-law model with a photon index T = 4.19^0 32 
(fis-so kcv = 2.4 x 10" 9 erg cm" 2 s~\ x 2 /dof=22/20). 

We extracted strictly simultaneous spectra from the XRT and 
BAT event lists (133-783 s since the BAT trigger; see Fig. [5) and 
we fit them in the 0.5-10 keV and 15-60 keV energy bands for XRT 
and BAT, respectively, by adopting the same models as those used 
for XTE J1739— 302. Factors were included in the fitting to allow 
for normalisation uncertainties between the two instruments, con- 
strained within their usual ranges (0.9-1.1). 

As is the case of XTE J1739— 302, for IGR J17544— 2619 a 
simple absorbed power-law model is inadequate to fit the broad- 
band spectrum (x 2 /dof=207 1/299). A significant improvement is 
obtained when considering the CPL and HCP models (Table[3](. 

The main difference in the x 2 -values between the CPL and 
HCP models is due to the residuals of the BAT spectrum above 
30 keV, which show in the former case a sinusoidal-like feature 
in the 15-30 keV region (see Fig. [6}. The inclusion of an absorp- 
tion edge in the CPL model improves the statistical result from 
X 2 /dof=383/298 to x 2 /dof=344/296, with best-fit parameters re- 
ported in Table [3] The F-test for discriminating among two differ- 
ent models (namely without and with the absorption edge) provides 
however a probability of chance improvement (PCI) of about only 
20%. The fact that these residuals are not observed in the HCP 
model besides on fact that in this case the spectrum is less smooth 
due to the presence of the E c term. The well constrained determi- 
nation of the latter parameter is actually indicative for the presence 
of a change of curvature in the X-ray spectrum around 15-20 keV. 

The models WABSx COMPTT yields x 2 /dof=339/297, and 
the inclusion of the absorption edge improves the fit to 
X 2 /dof=323/295, with however an even lower PCI (~ 35%) with 
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Figure 6. Residuals to data in units of a of the four models adopted to fit 
the continuum of IGR J17544-2619 (see Sect. [J/2]for details). 



respect to the case of a CPL, and best-fit parameters reported in 
Table H] 

Similarly to XTE J1739— 302, a satisfactory fit can be ob- 
tained by the sum of two blackbody spectra (see Table |4), and in 
this case the inclusion of the absorption edge only marginally im- 
proves the fit by A^ 2 ~ 7. 

Performing the spectral analysis with the COMPMAG model, 
we followed the same procedure of XTE J1739— 302. First, we 
considered the case of free-fall velocity profile fixing = 0.05, 
accretion column radius ro = 0.25 and NS albedo A = 1. 
In this case we obtain x 2 /dof=343/297, with an improvement 
Ax 2 ~ 20 when including the absorption edge at 20 keV. Also for 
IGR J 17544— 2619, the fit results to be insensitive to ro, with the 
X 2 -value substantially unchanged for ro increasing from 0.25 to 1. 

In the second case, namely = 0.2, the inclusion of the ab- 
sorption edge reduces similarly the fit by Ax 2 ~ 20, but we ob- 
served a decreases from x 2 /dof=363/295 to x 2 /dof=332/295 by 
setting ro = 0.25 and ro = 1, respectively. For higher values 
of ro, no further improvement of the \ 2 was observed, so that we 
consider it as our actual best-fit parameter. 

Unlike the case of XTE J1739— 302, the increased value of 
leads to a more significant decrease of the electron temperature, 
with kT c — 1.2+Q* keV, while the optical depth to first approxi- 
mation remains constant within errors, with r = 1.44j^Q'Jg. 

The reason for this is that a higher value of is compensated 
by the simultaneous increase of the accretion column radius ro 
from 0.25 to 1. 

In Fig. [6] we report as a summary, the residuals to the data 
in units of a of the four different adopted models above de- 
scribed, while Fig. [7] shows the absorption-corrected best-fit mod- 
els and residuals between the data and the model in units of a for 
BBODYRAD+BBODYRAD and EDGE X COMPMAG. 



4 DISCUSSION 

XTE J1739-302 and IGR J17544-2619 are considered the proto- 
types of the SFXT class and, as such, were systematically observed 



Energy (keV) 

Figure 7. Asborption-corrected EF(E) spectra, best-fit models and residuals 
between the data and the model in units of a for IGR J17544— 2619. Upper 
panel: BBODYRAD+BBODYRAD. Lower panel: EDGE X COMPMAG. 



with Swift, which caught several outbursts from both sources. The 
simplest way to parametrise the shape of the broad-band spectra 
is to fit them with phenomenological models typically used to de- 
scribe the X-ray emission from accreting pulsars in HMXBs, i.e., 
absorbed power laws, absorbed power-law models with a high en- 
ergy cut-offs and absorbed power-law models with exponential cut- 
offs. 

The advantage of this approach is that the fits yield easy-to- 
compare estimates of fluxes and luminosities for each outburst. 
Indeed, the data presented in this paper yield values of both high 
energy cut-off E c and e-folding energy E{ (see Sect. 13.11 and 
I3.2t > c onsistent with t he ones obtained prev i ously f or both sources 
( Sidolietal . 2009b a; Roma no et al.l201 ld l3, l20ich . The disadvan- 
tage is, clearly, that little physical insight can be obtained from 
such fits. In this paper we also apply the following more physically- 
motivated models: absorbed Comptonization models (COMPTT); a 
combination of two blackbodies with different temperatures and 
radii; and for the first time to SFXTs, the new COMPMAG model, 
recently developed by F12. 

There is an interesting aspect to discuss, which has not been 
yet fully faced in the accretion physics of SFXTs. If we focus on the 
results using the COMPTT and COMPMAG models, it can be shown 
that the electron density where Comptonization takes place is of the 
order of 10 19 cm -3 . Indeed, the Thomson optical depth is given by 



7 nig re, 



(1) 



where nig = n c /10 19 and re = i?/10 6 are dimensionless electron 
density and system length scale, respectively. If spectral formation 
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occurs close to the NS (re ~ 1-2), from the COMPTT best-fit values 
of r reported in Table[5]we obtain nig ~ 1. 

The possible presence of a strong magnetic field (B > 10 12 
G), with associated reduction of the Thomson cross-section <tt 
would actually require an even higher electron density, which 
somewhat compensates the lower cross-section value in order to 
maintain the same Comptonization parameter Y, and in turn the 
observed spectral index. 

This qualitative effect can be tested from the results of the 
COMPMAG model (Table [6j. The electron density for column ac- 
cretion, which is assumed in COMPMAG, is given by 



10 1 



mrg/3 z 



(2) 



where m = M / M^dd is the accretion rate in Eddington units, 
m — Mns/Mq is the NS mass in units of solar masses, ro = 
Ro/(R,Q W m) is the accretion column radius in units of the NS 
Schwarzschild radius, and /3 Z = V z /cis the accretion column ve- 
locity. We derived m from the best-fit value of r reported in Table|6] 
for IGR J17544-2619 and XTE J1739-302 (with r = 0.25 and 
= 0.05) and using equation (44) in F12 for both sources. Then, 
substituting /3 Z in equation ((2} with a value averaged over the ver- 
tical z-coordinate < /3 Z >, we obtain n c ~ 10 21 cm~ 3 in both 
sources, which is actually a factor about one hundred higher than 
that inferred from COMPTT. 

Then, we consider the mass flow rate across the magneto- 
sphere 



(3) 



where i? m is the magnetospheric radius, and Vi n is the infall veloc- 
ity at Rm, that is expected to be less than the free-fall velocity 



/2GM NS \ 1/2 
VB= {— —) ■ 

The accretion luminosity at the NS surface on the other hand is 

r „ GMnsM 



(4) 



Rns 



(5) 



Assuming a NS with mass M = 1.4M© and radius _Rns = 10 km, 
and combining equations $3$ and ((5J we obtain 

n a > 1.3 x lO 14 -^ cm" 3 , (6) 

rivg 

where L37 = L x /10 37 , rg = r m /10 8 , and 119 = iii n /10 9 , respec- 
tively. 

We consider a magnetospheric radius r m ~ 10 s cm, calcu- 
lated assuming the typical values of the magnetic field of an ac- 
creting pulsar in HM XBs and the typical stel lar wind parameters of 
OB supergiants jDavidson. & Ostrikeilll973l) . Additionally, for the 
infall velocity we assume for sake of simplicity V[ n « ve ■ 

With these prescriptions in mind and noting that both sources 
during the joint XRT/BAT observation have shown L x < 10 37 
erg s ~ 1 , from equation ((6) we find that at the magnetospheric radius 
n e > 10 14 cm" 3 . 

This value of the electron density at the magnetospheric radius 
is about seven orders of magnitude lower than that inferred from the 
best-fit parameters of COMPMAG as shown above, and thus cannot 
be simply a result of a simplified treatment of the problems such as 
pure spherical accretion assumed in equation ((6). 

However, there are at least to issues which are worth mention- 
ing. First, the electron density estimated from the best-fit parame- 
ters of COMPMAG is related to the region where the X-ray spectral 



formation is assumed to take place, which, in the case of COMP- 
MAG, correspond to cylindrical accretion column close to the NS 
surface with characteristic height-scale H ~ 1-2 /Zns- 

This scale-length is about 1/100 of the assumed magneto- 
spheric radius (« 100i? ns ), and for a density scaling as R~ 2 it 
would imply to first approximation a density increase from i? m 
.Rns about n e /n™ ~ 10 4 . Of course, matter channeling towards 
the magnetic poles may also play an important role in increasing 
the matter density in the region close to the NS surface. 

The second point to be considered is that if the NS orbital 
motion is supersonic, strong shock waves are expected to form ap- 
proximatively at the NS magnetosphere (bow shock region). The 
net result would be thus an increase of the electron density with 
respect to the value reported in equation (|6j. 

The Mach number of the NS is given by 

M w 10 2 Vs ns (~fZT c /fi)- 1/2 , (7) 

where V 8 nB is the NS orbital velocity in units of 10 s cm s _1 , 7 is 
the adiabatic index of the wind, /1 = mj /m p its molecular weight, 
Z its the charge state, and T c its temp erature in eV. For wind tem- 
peratures in the range 10 5 — 10 7 K jDucci et"al]|2009l ) , and ve- 
locities of the NS derived from the orbital parameters for both for 
XTE J1739-302 and IGR J17544-2619, M varies from a few to 
about 30, and if the shock is isothermal, p\j pi ~ M 2 . Actually, 
a density increase of a factor ~ 10 2 at the magnetospheric shock 
front, with respect to the value computed in equation ((6), would 
be sufficient to take into account the electron density derived from 
COMPMAG at the magnetic poles (see above). 

The physics of shocks in these systems is however rather com- 
plicated for the intrinsic three-dimensional nature of the problem. 
Indeed, if from one side the supersonic motion of the NS ensures 
the formation of a bow shock, on the other hand this discontinuity 
occurs in a region where the presence of a strong magnetic field 
plays an important role in determining the gas configuration. 

A pure dipolar magnetic field may lead to the formation of a 
bow shock which follows approximatively the shape of the mag- 
netic lines in the direction of motion of the NS, together with a 
channeling of matter towards the magnetic poles. This accretion 
geometry close to the NS could be approximated with a cylindri- 
cal accretion, suitable to be described with some accuracy by the 
COMPMAG model. However, if non-negligible multipole magnetic 
field components are presents, then matter can also (or mainly) be- 
ing accreted at the equator of the NS. 

As shown in the data analysis, the spectra of both 
XTE J1739-302 and IGR J17544— 2619 can be alternatively de- 
scribed by a two-component model consisting of a soft and hard 
BB (see Table [4]l . The temperature of the softer BB is consistent 
with that of the seed photons in both COMPTT and COMPMAG, 
while the temperature of the harder BB is comparable with the 
electron temperature of the above mentioned Comptonization mod- 
els. This means that the harder BB is actually playing the role of a 
Comptonization feature, but the quality of the high-energy data (see 
Figs.[4]and[7) does not allow to distinguish between an unsaturated 
(COMPTT and COMPMAG) or saturated (BB) Compton regime. 

The apparent radii of both BB components are consistent with 
emission regions of the order of a polar cap radius, thus pointing in 
favour of a very compact emission region. The implications for the 
accretion geometry in this context are that both the regions of the 
soft seed photons and of the Comptonized ones are visible. 

Summarizing, the current data presently do not allow us to dis- 
tinguish whether the 0.1-60 keV X-ray spectra of XTE J 1739 -302 
and IGR J 17544— 2619 form as a result of a single unsaturated 
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Comptonization process, or arise from two different and directly 
visible zones of seed and Comptonized photons. 

Us ing the calibration o f stellar parameters for Galactic O 
stars of lMartins et all d2005l) . we find Msg ~ 31 M , R SG ~ 

21 R Q for XTEJ1739-302, and Msg ~ 30 M , i? S G ~ 

22 R e for IGR J17544— 2619, respectively. If the orbital period 
of XTE J1739— 302 is 51.47 d, the minimum and maximum dis- 
tances of the NS from the giant companion are in the range 
d m in ~ 8 — 4.5 i?sG and d max ~ 9.8 — 13.4 i?sG for orbital 
eccentricity e = 0.1 — 0.5, respectively. These values change to 
rfmin ~ 3.2 — 1.8 Rsg and d max ~ 3.9 — 5.3 Rsa if the orbital 
period is 12.89 d (see Section[T). 

In the case of IGR J17544— 2619, on the other hand, the maxi- 
mum allowed eccentricity for its orbital period of 4.93 d is e = 0.4, 
and the ranges of minimum and maximum distances to the su- 
pergiant companion is d m i n ~ 1.6 — 1.1 -Rsg and d max ~ 
1.9-2.4 Rsg- 

The mass loss rates of the supergiant companions of 
XTEJ1739-302 and IGR J17544— 2619 have, within the uncer- 
tainties with which we know the masses, radii, effective tempera- 
tures and bolometric lumino sities of the two s tars, the same value 
of about 2 x 10~ 6 M Q yr _1 dVink et alj|2000h . 

Whatever the actual orbital period of XTE J1739 — 302 and the 
eccentricity of the two systems, the distance to the supergiant com- 
panion is lower for IGR Jl 7544— 2619. It is thus intriguing to ob- 
serve that the optical depth derived from the Comptonization mod- 
els is higher in the latter source than in XTE J 1739— 302 (see Tables 
[5]and|6]l, as one would expect for a closer orbiting system. 



5 CONCLUSIONS 

We have analyzed data from the outbursts of two sources, 
XTE J1739-302 and IGR J17544-2619, which are considered to 
be the prototypes of the SFXTs class. During the bright flare they 
reached peak luminosities of ~ 2 x 10 36 erg s _1 and ~5x 10 36 
erg s _1 , respectively. 

The presented Swift data do not allow us to discriminate 
whether the X-ray spectra are the result of a single unsaturated 
Comptonization process, or whether we observe both regions 
where seed photons are produced and Comptonization mostly takes 
place (two-BB model). 

The electron density in the region of the X-ray spectral forma- 
tion, computed to first approximation using the best-fit parameters 
of the COMPMAG model, is a factor about 1000 higher than ex- 
pected from the continuity equation at the magnetospheric radius. 
We propose that the formation of a bow shock at the NS magneto- 
sphere, due to its supersonic orbital motion around the supergiant 
companion, may increase the density enough to explain the differ- 
ence. This effect needs to be investigated in a future work by means 
of 3D magnetohydrodynamical simulations using the specifically 
devoted FLASH code. 

The possible feature around 20 keV which is observed using a 
one-component Comptonization model (see Fig. [7] bottom panel) 
is intriguing, as it can not be attributed to any particular shape of the 
BAT instrument effective area, however that fact that it is strongly 
dumped using a two-BB model and it has a low statistical signifi- 
cance do not allow us to draw firm conclusions. 

Nevertheless, this feature needs a further investigation with 
higher-resolution missions in the energy range around 20 keV, in 
particular with NuSTAR. If true, it would be the first measure ever 
obtained of the magnetic field intensity (~2x 10 12 G) in a SFXT 
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APPENDIX A: THE OPTICAL DEPTH RANGES OF APPLICABILITY OF THE COMPMAG MODEL 



The theoretical spectra of the COMPMAG model were reported by F12 for the case of optical depth r < 1. In this section we better explain 
how this impacts the spectral fitting analysis of real data. In particular, the fits of the broad-band spectrum of XTE J1739— 302 with the 
COMPMAG model (see Table |6) yield an optical depth formally less than 1, albeit within a factor of 3 from unity. It is worth noting that in 
COMPMAG this is the vertical optical depth of the accretion column and it is computed including an energy-independent correction to the 
Thomson cross-section, due to the presence of a magnetic field B > 10 12 G. More specifically, it is the opacity related to photons which 
diffuse across the field lines. The relation between the optical depth presented in Table|6]can be written in terms of Thomson optical depth as 
r « 1CT 3 tt, that in the case of XTE J1739-302 is r T ~ 300. 

After rearranging some terms, equation (35) in F12 can be written as 
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In the escape time prescription provided by BW07, the spatial diffusion of photons is described by 

le\\d 2 n /£^ 2 
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or, more clearly, by 
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where dry = dr — n e a\\dZ is the optical depth along the Z-axis and dr± — n e a±dr is the optical depth along the r-axis, perpendicular to 
the magnetic field, and A is the eigenvalue. The right hand side of equation lA2l can be approximated as 
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where c/tt = n^a^dZ and since the eigenvalue A 2 < 1, also the coefficient of n on the right hand side of equation dA3t should be less 
than 1, in order for diffusion approximation to hold. In particular, the coefficient Dy = ay / (3<7T| 2 ) takes into account the spatial diffusion of 
photons propagating along the magnetic field lines, while the coefficient T>± = l/(ron e aT±) accounts for the diffusion of photon travelling 
perpendicular to the field. 

Following the approach of B W07, from equation ((A3) it follows that the diffusion approximation for photons propagating in the direction 
of the magnetic field is valid if the corresponding Thomson optical depth tt > 10 . This can be actually the condition for a meaningful 
applicability of the COMPMAG model, together with the requirement that both the optical depths tu and t± need to be larger than unity. From 
equation jA3\ it is also evident that the ratio between the spatial diffusion in the Z-direction (2?y) and in the r-direction (V±) depends upon 
the choice of the column size (Zo/ro). 
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